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Abstract 

Many particle physics models attempt to explain the 130 GeV gamma-ray feature that the Fermi-LAT 
observes in the Galactic Center. Neutralino dark matter in non-minimal supersymmetric models, such as 
the NMSSM, is an especially well-motivated theoretical setup which can explain the line. We explore the 
possibility that regions of the NMSSM consistent with the 130 GeV line can also produce the observed baryon 
asymmetry of the universe via electroweak baryogenesis. We find that such regions can in fact accommodate 
a strongly first-order electroweak phase transition (due to the singlet contribution to the effective potential), 
while also avoiding a light stop and producing a Standard Model-like Higgs in the observed mass range. 
Simultaneously, CP-violation from a complex phase in the wino-higgsino sector can account for the observed 
baryon asymmetry through resonant sources at the electroweak phase transition, while satisfying current 
constraints from dark matter, collider, and electric dipole moment (EDM) experiments. This result is 
possible by virtue of a relatively light pseudoscalar Higgs sector with a small degree of mixing, which yields 
efficient s-channel resonant neutralino annihilation consistent with indirect detection constraints, and of the 
moderate values of /i required to obtain a bino-like LSP consistent with the line. The wino mass is essentially 
a free parameter which one can tune to satisfy electroweak baryogenesis. Thus, the NMSSM framework can 
potentially explain the origins of both baryonic and dark matter components in the Universe. The tightness 
of the constraints we impose on this scenario makes it extraordinarily predictive, and conclusively testable 
in the near future by modest improvements in EDM and dark matter search experiments. 
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I. INTRODUCTION 



In the search for signatures from the annihilation (or the decay) of dark matter particles, a 
gamma-ray line in the multi-GeV energy range has long been considered a Holy Grail. Given that, 
in the weakly interacting massive particle (WIMP) paradigm. Galactic dark matter is virtually 
at rest, the pair annihilation of two particles into a final state consisting of two photons would 
produce a monochromatic line with an energy exactly corresponding to the particle dark matter 
mass (or to half its mass in the case of decay). The advent of the Fermi gamma-ray Large Area 
Telescope (LAT) heralded promise of potentially delivering this smoking gun signal, which would 
then serve as a beacon for further searches to close in on a well-defined particle dark matter mass. 

Despite a null result presented by the LAT collaboration in Ref. independent scholars an- 
alyzed the Fermi data ernploying optimized signal-to-noise regions, unveiling a tantalizing excess 
localized around 130 Ge\|j and originating from regions including the Galactic center 3j,|j]. Subse- 
quent independent analyses confirmed the original claim, typically attributing an even larger level 
of confidence to the discovery of a monochromatic line in the Fermi-LAT data from the center of 
the Galaxy [s]. 

Understandably, the discovery of the line spurred a great deal of interest in the community: 
a feature in the Earth limb photon events at the same energy was found, albeit with a much 
lower statistical significance 5|; despite significant efforts in pinpointing possible instrumental or 
environmental effects that could explain the excess (see e.g. Ref. ^), at present the line feature 
appears statistically significant enough to deserve serious consideration. 

From a model-building and phenomenological standpoint, the 130 GeV line poses interesting 
challenges: with default choices for the dark matter density profile in the Galaxy, the required pair- 
annihilation cross section for dark matter (at rest, i.e. at "zero temperature") into two photons is 
about {av)^^ ~ 10^^'^cm'^/s, much larger than would be expected by suppressing by a factor 
the pair annihilation cross section expected for WIMP thermal production in the early universe. 
Even more problematic is the absence of a continuum gamma-ray signal accompanying the line in 
the region where the line is detected. This poses the question of how to suppress final states that 
would generously produce e.g. neutral pions from hadronization showers of strongly interacting 
particles, or inverse Compton or bremsstrahlung photons from charged leptons. 

Simple paradigms for WIMP dark matter fail at explaining the needed features of the 130 
GeV line. For example, neutralinos within the minimal super symmetric extension of the Standard 



^ Recent re-analyses with reprocessed data using "Pass 7 Clean" events put the line at 135 GeV [3, but nothing 
qualitative changes in the present discussion, where we will assume the line is at 130 GeV. 
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FIG. 1. The dominant diagram leading to the two-photon pair- annihilation of neutralinos in the NMSSM 
scenario under consideration in this study. 



Model (MSSM) feature large suppressions in the pair annihilation into two photons with respect to 
any other final state, and the required large rate for neutralino pair-annihilation into two photons 
cannot be accommodated with the right thermal relic abundance Q]- 

A simple extension to the field content of the MSSM, however, allows for an interesting caveat 
to both shortcomings mentioned above, as first realized in Ref. [s!]: within the next-to-MSSM (or 
NMSSM, hereafter), an s-channel resonant contribution exists to the annihilation cross section 
arising from the diagram shown in Fig. [H where two approximately 130 GeV bino-like neutralinos 
annihilate into a singlet-like pseudoscalar Ai, which then decays into photons via a chargino loop. 
For rriAi ~ 260 GeV, the process is resonant and the resulting cross-section can easily satisfy 
(cTw)^^ ~ 10^^ ''cm^/s as required to produce the observed line [31]. 

The NMSSM possesses the interesting additional possibility of naturally realizing a mechanism 
known as electroweak baryogenesis to produce the observed baryon asymmetry of the universe 
(BAU) at the electroweak phase transition (EWPT) (for a recent review, see Ref. 9|]). The NMSSM 
framework, in fact, accommodates tree-level cubic couplings entering the relevant scalar effective 
potential driving the EWPT needed to produce a sufficiently strongly first-order phase transition 



(this is in turn needed to prevent wash-out of the generated baryon asymmetry in regions of bro^ 



Ken 



electroweak phase), as realized a long time ago 10|, lUl] and reinforced in recent analyses [12 1 
(see Refs. [l3 . [l^ for similar arguments in related models). Additionally, the NMSSM, like the 
MSSM, possesses enough room to host the level of CP violation needed for baryogenesis while 
being consistent with constraints from the non-observation of electric dipole moments (EDMs). 
In the present study, we argue that the NMSSM can simultaneously accommodate: 



1. a thermal dark matter candidate that can produce the 130 GeV line while being consistent 
with constraints from other gamma-ray observations and direct detection searches; 
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2. a Higgs sector consistent with the recent LHC findings 




3. a strongly first-order phase transition as needed by electroweak baryogenesis (for which we 
calculate in detail the effective finite temperature potential); 

4. the generation of the observed baryon asymmetry of the universe at the EWPT, while being 
consistent with constraints from EDMs. 

Requiring all four conditions above forces us to very special corners of the theory's parameter 
space: the goal of our study is not to explore exhaustively the NMSSM parameter space but, rather, 
to outline the general implications for the theory parameter space of the four requirements above, 
and to draw predictions from the regions of parameter space that do satisfy these requirements. As 
a result, we do not concern ourselves with issues of fine-tuning but, rather, we produce a detailed 
set of predictions that put this framework for the origin of baryonic and dark matter on very 
testable grounds. At the same time, we provide benchmarks for corners of the NMSSM theory 
parameter space where all conditions listed above may be fulfilled. 

This paper is organized as follows: in Sec.[TI]we outline the NMSSM parameter space, detail the 
neutralino and Higgs sectors, and discuss the phenomenological constraints we implement; Sec. IIIII 
discusses the nature of the electroweak phase transition and the constraints that a strongly first- 
order transition places upon the parameter space; in Sec. |IV] we discuss the computation of the 
baryon asymmetry; we conclude in Sec. |V] 



II. A 130 GEV LINE IN THE NMSSM 



To begin, we review the NMSSM setup, and show how it is possible to hone in on parameters 
consistent with the 130 GeV gamma-ray signal and with a broad set of additional phenomenological 
constraints. We follow closely the strategy outlined in Refs. P, and consider the simplest 
incarnation of the NMSSM with a scale- invariant, Za-symmetric superpotential: 

W = WMSSM\^.=o + XSHuHd+'^S^ (1) 

where hatted quantities denote the corresponding superfields, and where 5 is a gauge singlet. The 
soft super symmetry-breaking Lagrangian is given by 

- £'"f' = -CZssM + \S\^ + (^XAxSH^Hd + ^kA^S^^ + h.c. (2) 
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After electroweak symmetry breaking (EWSB), the Higgs and singlet fields obtain vacuum expec- 
tation values (vevs) of {Hy) = Vu, {Hd) = Vd, and (S) = Vg- As in the MSSM, we denote the 
ratio of the SU{2) Higgs vevs as tan/3 = Vu/vd- The singlet vev generates an effective /i-term in 
the superpotential given by ^ = Xvg- We assume that A,Us G M so that ^ is real and there is no 
CP-violation at tree level in the Higgs sector. While CP-violating effects can enter at one-loop 
from gaugino interactions if we allow Mi^2 to carry a complex phase, we neglect these contribu- 
tions when considering radiative corrections to the Higgs sector, since these effects are typically 
sub-dominant. The six parameters A, k, Ax^ A^, ^ and tan/3 then determine the tree-level Higgs 
spectrum after minimizing the scalar potential and solving for the SUSY-breaking Higgs masses. 

At this level, deviations from the spectrum of the MSSM originate from the singlet superfield 
in the superpotential, and are crucial in order to obtain a neutralino consistent with the 130 GeV 
gamma-ray signal (without an associated continuum gamma-ray background), with a 125 GeV 
Higgs, and with successful electroweak baryogenesis. Specifically, the present set-up contains one 
each of additional neutral CP-even and CP-odd states which enter into the respective Higgs mixing 
matrices. Complete expressions for the various relevant mass matrices in the NMSSM which match 
our conventions can be found in, e.g., Ref. [3]. 

The pseudoscalar mass matrix will be of particular importance; its elements are given, to one- 
loop order, by [3] 

A. 2 X , Aan/3(Q) cot/3(Q)\ 

-^1,11 = (^A + ^vs) y-^^ + ) 

M%22 = ^X^iVu{Q)vd{Q) + A^A^^^i^^^^^ - SkA^Vs (3) 

■Mp,12 = A ^ 2^ J (^A - ^KVs) , 

where Q is the relevant SUSY energy scale; Vu,d{Q) and tan/3((5) are the Higgs vevs and tan/3 
at the scale Q; and Zh^^{Q) are wave-function renormalization factors. The matrix Mp can be 
diagonalized to obtain the pseudoscalar mass eigenstates Ai and A2. As we discuss below, in the 
present setup Ai must be singlet-like; the state A2 will therefore correspond to an MSSM-like 
pseudoscalar Higgs boson. 

In addition to the new degrees of freedom in the Higgs sector, there is an additional Weyl 
fermion (the "singlino", S), corresponding to the fermionic component of the singlet superfield S. 
This fermionic degree of freedom enters into the neutralino mixing matrix, whose components are 



given at tree level by 
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Here, we shall consider the case in which the baryon asymmetry is sourced by CP-violation in the 
higgsino-gaugino sector ^. The masses in Eq. are therefore generically complex- valued. We 
will further restrict ourselves to the case of a single complex physical phase, in the wino mass M2, 
with all other parameters reaj^. This results in CP-conservation at tree-level in the Higgs sector. 
Since in our construction the LSP is bino-like throughout all of the parameter space we consider, a 
CP-violating phase in Mi would produce large effects on the calculation of the various dark matter 
properties; we therefore impose Mi G M. Eq. Q is diagonalized by the unitary complex matrix A^: 



(5) 



and the neutralino masses are given by 

diag (m^o, mjo, m^o, m^o, mj.,) = M'IM^o- (6) 

The five neutralinos are admixtures of i?, W , Hu^d, and S, the lightest of which will be the lightest 
supersymmetric particle (LSP) in our setup. The chargino mass matrix is simply that of the 
MSSM, again with a possible complex phase in the wino mass entry, yielding the mass eigenstates 

Motivated by the lack of a SUSY particle discovery at the LHC, we will assume that all sfermions 
are heav}0, with rrigj > 1.5 TeV. This effectively decouples them from any processe of interest here. 
As a result, to determine the properties of neutralino dark matter, the electroweak phase transition, 
and the CP-violating sources for electroweak baryogenesis in the present set-up, one must specify 
the following nine NMSSM parameters: 



A, K, Ax, ^, tan/3. Mi, IM2I 



Note that the physical phase we consider here effectively corresponds we to the phase 
Ref. Q 



= arg(M2). (7) 
; arg(/iM2 6*), see e.g. 



Note that the authors of Ref. [3| considered rather light sleptons to account for the possible discrepancy of the 
muon .g — 2 with the value predicted by the SM. However, in the present case, such light sleptons can result in 
large one-loop contributions to the electric dipole moments inconsistent with the constraints discussed in Sec. II V Cl 
barring cancellations. 
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As we argue below, many of these parameters are tightly constrained by the phenomenological and 
observational constraints we impose, in particular by requiring a 130 GeV gamma ray line from 
resonant neutralino annihilation consistent with other particle and dark matter searches. 

Throughout this study, we will assume that the large required pair-annihilation cross-section 
into two photons, (av)^^ > 10^^'' cm^/s, arises from the on-resonance s-channel annihilation of 
neutralinos into Ai, which in turn couples to two photons through a chargino loop (see Fig.[T|). The 
dominant contribution to the thermally averaged cross-section for this process at zero temperature 
is given by 



(av) 



77 



167r3 



E 



i=i,2 'i'fU'Zo \ 4m^o 



m. 



where the function F{a,b) is defined by 



F{a,b) 



1 dx 



log 



Aax^ — 4ax + b 



(8) 



(9) 



and the couplings Qa o , g ^ ± depend on the neutralino, chargino, and CP-odd Higg s diagonalizing 

Iq . appropriately 



matrices. To compute these couplings, we use the Feynman rules found in Ref. 
modified to match our conventions for the neutralino and chargino matrices, which contain complex 
mass entries. This cross-section is plotted as a function of rriAj^ for a particular choice of parameters, 
in Fig. O which clearly shows the narrow resonant structure. 



A. Suitable Higgs and Neutralino Sectors 



Given our set-up, we can elucidate the parameter space regions capable of producing the gamma- 
ray line while satisfying all other dark matter and particle physics constraints. As we show below, 
requiring a 130 GeV line from resonant neutralino annihilation restricts the NMSSM parameter 
space to a narrow region in which we can study electroweak baryogenesis and the electroweak phase 
transition, in addition to producing unambiguous predictions for several experimentally observable 
quantities, such as electric dipole moments and dark matter detection rates. 

In general, the properties associated with the neutralino LSP depend sensitively on the details 



of the various parameters invo^ 



points discussed in Refs. 



ved; this can be appreciated by considering the different benchmark 



17l |. For example, the annihilation cross-section into photons, Eq. 



is strongly affected by the mass splitting 



2m^o 



as shown in Fig. [2j Correspondingly, other 



resonant processes, such as the s-channel neutralino pair annihilation into bb through ^i, also 
depend on the mass difference. The details of the various resonant channels significantly affect 




FIG. 2. The zero-temperature thermally-averaged cross-section times velocity for neutralino annihilation 
into two photons as a function of the singlet-like pseudoscalar mass for the EWPT benchmark point 
discussed m Sec. \TE- X = 0.75, k = 0.45, tan/3 = 1.7, Ax = 545 GeV, = -88 GeV, n = 275.8 GeV, 
Ml — 143.5 GeV, and AI2 = 635.5 GeV. The red dashed line indicates the lower hound on (crw)^^ required 
to produce the 130 GeV Fermi line. Note that decreasing Mi (thereby increasing fi) will narrow down the 
resonance. 



both the zero-temperature and the finite-temperature annihilation cross sections (the latter being 
relevant for the calculation of the thermal relic density of dark matter) . The amplitudes associated 
with these processes can however be tuned so that the neutralinos produce a 130 GeV gamma-ray 
line while satisfying all other indirect detection and relic density constraints, as we show here. 

Since we will be concerned with properties of the electroweak phase transition and baryogenesis 
which do not depend sensitively on the details of the resonance, it is sufficient, for our purposes, 
to consider the simple parameter choice Ai = 2m^o = 260 GeV and proceed to consider the 
implications for electroweak baryogenesis (a slightly off-resonance value would not at all affect 
the electroweak phase transition or the resulting baryon asymmetry). From this starting point, 
we shall dial in the various parameters point-by-point to satisfy all of the phenomenological and 
observational constraints we describe below. 

First and foremost, besides requiring the desired neutralino annihilation structure, demanding 
a 130 GeV LSP neutralino and the associated 260 GeV singlet-like Ai, we require a 125 GeV SM- 
like Higgs, in accordance with recent experimental findings from the LHC collaborations [isl. 
Given our parameter space, the requirements on the bino-like LSP and on Ai lead us to vary Mi 
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and Ax in the range 

135GeV < Ml < 145GeV 

(10) 

150GeV < < eOOGeV. 

For each point in the Mi, Ax parameter space, we use the fohowing strategy to choose values for 
the seven remaining parameters: 

1. To obtain a Higgs mass of 125 GeV in the NMSSM without excessive tuning in the stop 
sector requires relatively large A and small tan /3, as seen from the tree-level inequality: 

2 ^ f , 2A^sin^2/3 ^ 2 

m^^^ < cos 2(3 + — o 2 ^z- (11) 

V Si ~r 92 / 

We take tan (3 in the range 1.7 < tan /3 < 1.8. In principle A can be either positive or negative. 
We focus on positive A and consider 0.6 < A < 0.8 (see, e.g. Ref. for a discussion of the 
case of A < 0). For |A| much smaller than this value, one must rely heavily on stop loops 
to raise the Higgs mass. Also, A determines the coupling of neutralinos to Ai, as well as 
the coupling of Ai to photons, and so for much smaller |A| the neutralino annihilation cross- 
section into photons is suppressed. For values A > 0.7, A becomes non-perturbative below 
the GUT scale; this can be remedied by including higher-dimension operators resulting from 
integrating out new physics which enters below the GUT scale. 

2. The pseudoscalar Ai must be predominantly singlet-like to be compatible with indirect 
detection results. The amount of mixing between Ai and the MSSM-like CP-odd Higgs A2 
is governed by Ai p,i2 in Eq. ([3]) and is minimized for 

(12) 

Given the relatively large values of A we consider, we take k > 0.3. For a given choice of 
K, the Ai — A2 mixing will vary point-by-point in the parameter space under consideration. 
Therefore in some regions of parameter space the lightest pseudoscalar can obtain a large 
branching ratio into fermions and be incompatible with indirect detection constraints for a 



given mass difference 



niA, — 2m^o 



. As mentioned above (and discussed in more detail in 
Sec. IIIBp . one can typically dial in the details of the resonance to satisfy these constraints 
for a given point, however the BAU does not depend sensitively on this tuning. 



■* We will in fact assume that this is the case for our benchmark EWPT point which features A = 0.75. 
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3. To obtain a lightest neutralino mass of 130 GeV, we must fix /i and M2 or, equivalently, ^ 
and A appropriately, where we define the quantity A via 



When considering CP-violation in Sec. IIVI we will typically set the CP-violating phase (p to 
its maximal value, svcKp = 1, in our calculations to show the maximum extent of the EWB 
parameter space, although viable regions will typically have phases of O(10~^). In calculating 
the baryon asymmetry, A will govern the strength of the resonant CP-violating source. In 
considering the higgsino-gaugino CP-violating sources we will typically take A = as an 
optimistic EWB scenario. Given a particular choice of A and 0, we fix by diagonalizing 
Eq. dH) and solving for ^ such that m^o = 130 GeV (note that we can rewrite Vg = /^/A). 
This procedure fixes all the relevant parameters in the neutralino and chargino sectors. 

4. Finally, to obtain a large photon annihilation cross-section, we need the annihilation channel 
XiXi ~^ ^1 to be near resonance at T = 0, which implies m^i ~ 260 GeV. As discussed above 
and shown in Fig. [2l there is a narrow (< 1 GeV) window for which {crv)^^ is large enough 
to be compatible with the line. Since the properties of the electroweak phase transition 
and baryogenesis are not sensitive to the precise value of rriAi, choose to sit exactly on 
top of the resonance, i.e. enforce mA_^ = 260 GeV, by diagonalizing Eq. ([3]) and solving 
for the appropriate value of Af^. Therefore, at each point in the parameter space, (cr^^v) > 
10~^^cm^/s. Once again, the precise mass splitting between Ai and the LSP can typically 
be tuned point-by-point to produce the line while providing the correct relic density and 
satisfying the other indirect detection constraints as described below. 

The strategy outlined above is useful to automatically select the regions in the NMSSM producing 
the tentatively observed 130 GeV gamma-ray line, and provides an efficient way to study the 
properties of electroweak baryogenesis in these regions by exploring the remainder of the parameter 
space. Note that we are not concerned with tuning or naturalness in this scenario, since we have 
narrowed in on this region by demanding consistency with the (tentative!) observation of a gamma- 
ray line which we postulate to be associated with dark matter pair annihilation. 

We shall now use our suitably selected Higgs and neutralino sectors to close in onto electroweak 
baryogenesis in regions of the NMSSM producing a 130 GeV line. However, we first comment 
further on the impact of various other dark matter and particle physics constraints on the parameter 
space under consideration. 




(13) 
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B. Phenomenological Constraints 



The NMSSM parameter space of interest features relatively light neutralino, chargino, and 
Higgs sectors and is thus quite constrained on multiple fronts. Here we highlight the most im- 
portant constraints on the parameter space and consider their impact on our current set-up. We 
use NMSSMTools 3.2. iQ and MicrOmegas 2.4.5[2l| to calculate the various cross-sections and 
quantities of interest. We summarize in Fig. [3] the impact of the constraints we consider here (and 
that we discuss in detail below) on the relevant parameter space, for the particular choice A = 0.6, 
K = 0.32, and tan /3 = 1.8 as an illustrative example. In these calculations, we take M2 to be 
real; since the LSP has only a very small wino component across the parameter space, and since 
the other neutralinos and charginos are significantly heavier than the lightest neutralino, the DM 
constraints will be largely unaffected by allowing M2 to be complex. The Higgs couplings are also 
insensitive to 6. 



1. Indirect Dark Matter Detection and Thermal Relic Density 

Indirect detection places important constraints on the parameter space in question. In con- 
sidering raA^ ~ 2m^o, there will also be a resonant tree-level neutralino annihilation channel into 
quark- antiquark, and especially 66, final states, eventually leading to gamma rays via hadroniza- 
tion producing neutral ions. The lack of an excess of gamma-rays associated with this emission 



puts constraints on the branching ratio for neutralino pair-annihilation into, e.g., 66 



23] . As men- 



to obtain both 



221 of the 



tioned above, however, one can generally dial in the mass splitting 
{(T^^v) > 10"^^ cm'^/s and {(Jfjiv) < 10"^^ cia^/s as required by Fermi observations 
diffuse gamma ray background (see e.g. the benchmark point in Table [J). Additionally, neutralino 
annihilation into W~^W~ will receive a contribution at tree- level from the pseudoscalar channel; 
however, this contribution also typically falls well beneath the 10"^^ cm^/s bound from Fermi 
by adjusting m^^. Consequently, this tuning allows one to satisfy all continuum gamma-ray con- 



straints 



23] while reproducing the observed intensity of the 130 GeV line, something that cannot 



be done in the MSSM. The parameter space we consider for electroweak baryogenesis can thus 
be dialed in to agree with indirect detection results without drastically affecting the details of the 
electroweak phase transition or the generation of the baryon asymmetry. 

Similar reasoning applies to the DM thermal relic abundance. For Xi to be a suitable thermally- 
produced dark matter candidate, it must be compatible with the bounds on the relic density from 
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FIG. 3. An example of the NMSSM parameter space for successful electroweak baryogenesis and a 130 GeV 
gamma-ray line. Here we take A — 0.6, k, — 0.32, tan/3 — 1.8 and A = (so that the CP-violating sources 
are on resonance) , while the rest of the parameters are chosen as described in Sec. \II A\ to be consistent 
with the Fermi line. The gray shaded region is excluded by the XENONIOO 225 live day results, calculated 
with the default settings in MicrOmegas. Red shaded regions are excluded by measurements of the Higgs 
mass (although these regions can be shifted around by changing e.g. the squark masses). The orange shaded 
region is excluded by the non-observation of an electric dipole moment of the electron. The blue contours 
correspond to points consistent with the observed baryon-to- entropy ratio of the universe for different values 
of the CP-violating phase 4>. 



WMAP7 2S['- ^DKih^ = 0.112 lb .011. While at zero-temperature the neutrahno sits very close to 
the pseudoscalar resonance, at the freeze-out temperature Tf.o. ~ m^o /20 w 6.5 GeV, the resonance 
is shifted higher by about 10 GeV for the case of m^i^ = 260 GeV. This can be seen by evaluating 
the thermally-averaged center-of-mass (C.O.M.) energy, (s), at T = Tf ^., given by 



(s) ~ 4mJo + 6m^oTf.o. ~ 270 GeV. 



(14) 



However, in evaluating (erf) at Tf.o., one integrates over center-of-mass energies, and hence effec- 
tively picks up contributions from the resonances, which decrease as one moves (s) further away 
from 4m^(, . Therefore, as is the case for the zero-temperature cross-sections, by dialing in the de- 
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tailed neutralino and pseudoscalar masses, as well as the Ai — A2 mixing, one can typically achieve 
a total annihilation thermally averaged cross-section of (av) ~ 3 x 10^^^ cm^/s required to obtain 
the correct relic density^ 
Previous studies 



17| have relied on a sizable higgsino component in the LSP to drive the 
relic density down. However, this requires small values of /i which are difficult to reconcile with 
the most recent direct detection constraints, except in the case of cancellations which can occur 



for negative ^ as exploited in Ref. [13] (we have found it difficult to achieve a strongly first-order 



36 possible). Another 
with mass near 130 



EWPT consistent with the 130 GeV line for the < case, but it may still 
possibility is to open a co- annihilation channel by e.g. allowing a light stau 
GeV to drive the relic density down. Light staus are not yet significantly constrained by LHC 
searches and, interestingly, they could provide an explanation of the enhanced Higgs diphoton rate 
as observed by ATLAS, albeit for large tan/3 (see e.g. Ref. [25|). We do not pursue these avenues 



further, but emphasize that we find that the relic density (and the zero-temperature neutralino 
annihilation cross-sections) can be made to agree with observations in this scenario by tuning or 
other mechanisms that do not significantly affect the properties of the EWPT nor the calculation 
of the baryon asymmetry. Consequently, we do not focus on the detailed bounds from indirect 
detection or the thermal relic abundance point-by-point in our present study of EWB in this 
scenario, but we do emphasize that these constraints can all be met in principle, as illustrated by 
a worked-out example in the EWPT benchmark point we show explicitly in Table HI 

2. Direct Detection 

Unlike the case of indirect detection and relic density constraints, the bounds from DM direct 
detection (i.e. the scattering of the lightest neutralino off of nucleons) do not depend sensitively 
on the details of the resonance, but rather on the composition of the lightest neutralino. This in 
turn depends on Mi: larger values of Mi require smaller values of /j, to obtain m^o = 130 GeV and 
consequently enhance the spin-independent neutralino-proton cross section. 

We require that the LSP satisfy the current upper bound from XENONIOO for a 130 GeV 
WIMP for the spin-independent cross-sectioro, asi ^ 3 x 10^^ pb [26]. We show the impact 
of this constraint on our parameter space in Fig. O points excluded by XENONIOO are shown 
in the gray shaded region. These bounds are computed assuming default values for the various 



^ Of course with CP-violation in the gaugino sector one must verify that such a hght slepton satisfies constraints 
from EDMs. 

^ We also consider the bound on the spin-dependent cross-section, but the corresponding constraints are much 
weaker than those on asi in our scenario 

13 



underlying parameters, such as the quark content of the nucleon, local distribution of dark matter, 
etc. We employ the MicrOmegas 2.4.5 



21[ package for the calculation of the relevant scattering 
cross section, and employ the default parameters thereof. As expected, points with smaller fi 
values, and hence a larger higgsino component in ^-re ruled out. 

We note here that the exclusions are somewhat stronger than those reported in Ref. jsl due to 
the release of the 2012 XENON results (and consequently the window for is somewhat more 
constrained than that in Ref. [8]). Since these limits depend on parameters affected by significant 
uncertainty, they should also be taken with a grain of salt. For example, by considering the strange 
quark content of the nucleons near the end of the error bars from Ref. 27]] {cTt^n = 39 MeV, ctq = 43 
MeV), one can push the XENON limits out to allow Afi up to ~ 145 GeV consistent with the 2012 
XENONIOO results (see e.g. the EWPT benchmark point in Table IH). 



3. Higgs Constraints 

The lightest CP-even Higgs in our scenario is SM-like. We require that 124 GeV < m/^^ < 127 
GeV, in agreement with results from ATLAS [3] and CMS |l5|. The region of parameter space 
incompatible with these results is shown in Fig. [3] by points within the red shaded regions. We have 
also checked against constraints from hi — t- bb, tt, etc. as implemented in NMSSMTools 3.2. ij^]. 
The couplings of hi to the various SM fermions and gauge bosons all fall within ~ 3% of the 
corresponding SM predictions, hence well within experimental limits. 

The lightest CP-odd Higgs must also be compatible with collider searches. In particular, we 
verified that the couplings of Ai to bb, tt are small compared to that of the SM-like Higgs for 
compatibility with LHC results. In the parameter space under consideration, we find that the 
couplings of Ai are at most of order 1% of the SM Higgs couplings. 



4-. Other Considerations 

There are several other constraints which are in fact satisfied over nearly all of the parame- 
ter space we consider. Constraints from LEP on light charginos are everywhere satisfied, since 
charginos are always heavier than the 130 GeV LSP. Also, constraints from i?-physics, as imple- 
mented in NMSSMTools 3.2.1, do not constrain the parameter space since we consider small values 
of tan/3. Finally, we have also verified the absence of unphysical global minima of the effective 
potential for all points we consider, as well as the absence of Landau poles below the GUT scale. 
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with the exception of the EWPT benchmark point, for which we take A = 0.75. As discussed 
above, this issue can be remedied with the modest assumption of new physics entering below the 
GUT scale. 

In summary, Fig. [3] shows that there exist regions of NMSSM parameter space consistent with 
a 130 GeV gamma-ray line, a 125 GeV SM-like Higgs, and which can satisfy all relevant dark 
matter and experimental particle physics constraints. We can now proceed to investigate the 
phenomenology and properties of electroweak baryogenesis in these regions. 



III. THE ELECTROWEAK PHASE TRANSITION 



Successful electroweak baryogenesis requires a strongly first-order electroweak phase transition. 
In the absence of a strongly first-order transition, SU{2) sphaleron processes, which provide the 
necessary baryon number violation, are unsuppressed in the broken electroweak phase and tend 
to wash out any existing generated baryon asymmetry. The strength of the phase transition can 
be parametrized by the order parameter (p{Tc)/Tc, where Tc is the critical temperature, defined 
as the temperature for which the symmetric and broken phases are degenerat^j- To prevent 
sphaleron washout requires ip{Tc)/Tc > 1, which we take as the definition of a "strongly first- 
order" transitio As we will show in this section, this requirement can be readily satisfied in the 
region of the NMSSM compatible with the 130 GeV gamma-r ay li ne and without relying on a light 
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stop squark, as is instead typically required in the MSSM 

The strength of the electroweak phase transition is governed by the finite-temperature effective 
potential, which comprises several parts: the tree-level scalar potential, zero-temperature quan- 
tum corrections, finite-temperature quantum corrections, and thermal mass terms. The tree-level 
potential comes directly from the superpotential (Eq. ([T])) and the soft supersymmetry-breaking 
terms (Eq. 

VoiK, hd, s) = ^(5? + 92) {hi - hlf + ^K^s^ _ \\^,s^huhd + ^A^ [hlhl + [hi + hi)) 

, Q \ A 7 7 ^ 111 ^ 111 ^ 22 /-IP-\ 

+ —i^A^s - —XAxshuhd + -m^h^ + -m^h^ + -m^s . (15) 
The fields hu, h^, and s are defined by 





/f. = ^| °|;/f. = -^l J;S=-^. (16) 



^ Note that this quantity is not gauge invariant, see e.g. the discussion in Ref. [28l. |29 
* More precisely, one should actuaUy consider the system at the nucleation temperature, T„. However, the amount 

of supercooling in this model is small, and for simplicity we assume that r„ ~ Tc as in previous work. 
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We assume that the scalar fields are real at all temperatures, and we do not consider charged vacua 
(although we do ensure that the potential is stable in the charged and imaginary directions). 



Using MS renormalization, the one-loop zero-temperature quantum corrections are 

(17) 



where are the (possibly negative) field-dependent mass-squared values, Ui are their associated 
number of degrees of freedom, A is the renormalization scale, and c = ^ for the transverse polar- 
izations of gauge bosons while c = | for their longitudinal polarizations and for all other particles. 
The plus and minus signs are for bosons and fermions, respectively. The sum over the relevant 
particles i include all standard model particles (although we ignore fermions lighter than the bot- 
tom quark), the physical Higgs and other scalar particles, their associated Goldstone bosons, the 
neutralinos and the charginos. We work in Landau gauge where the ghost bosons decouple and 
need not be included in the spectrum. The one-loop potential contains explicit gauge-dependence 
which cancels with the implicit gauge-dependence of the vevs at every order in h (for recent dis- 



28|, 



2g 



331]). As is common 



cussions of gauge dependence in effective potentials, see e.g. Refs. 
practice, we do not consider the effects of the implicit gauge-dependence, and therefore our results 
will contain gauge artifacts. However, our primary purpose in examining the effective potential is 
to estimate whether or not a first-order phase transition is possible, and for this purpose a rough 
calculation with gauge-dependence is acceptable. 

We calculate the neutralino masses from Eq. @ above. The scalar mass matrix is given by taking 
the second derivative of the tree-level potential, but including CP-odd and charged directions. 
This yields a block-diagonal 10 x 10 matrix, with blocks consisting of CP-even states (3 degrees of 
freedom), CP-odd states (3 degrees of freedom), and two blocks of charged Higgses (4 degrees of 
freedom) (see Appendix [A] for details). 

The finite-temperature contributions are 

Fi(r>0) = V^{T = 0) + ^Y.n.J± , (18) 

where 

J±(x2) ^^j^dy y2 log (l ^ e-^/^) (19) 

and again the upper (lower) signs correspond to bosons (fermions). At high temperature, the 
validity of the perturbative expansion of the effective potential breaks down. Quadratically di- 
vergent contributions from non-zero Matsubara modes must be re-summed through inclusion of 
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thermal masses in the one-loop propagators SJ, l35| . This amounts to adding thermal masses to 



the longitudinal gauge boson degrees of freedom and to all of the scalars (see Appendix lAj) . 
The full one-loop effective potential is 

V{K,hd,s,T) = Vo{hu,hd,s) + Vi{T = 0) + ^^niJ± (^^^ 



(20) 



where the masses mf are field-dependent and include thermal mass corrections. 

The important qualitative feature of the finite-temperature contribution is that it lowers the 
effective potential anywhere mj/T'^ is small. To get a strongly first-order phase transition, we need 
to sharply lower the potential near the symmetric phase without significantly lowering it in the 
broken phase so that the two phases may be degenerate with a sizable barrier. Therefore, a strongly 
first-order transition demands either numerous heavy field-dependent particles (such that they are 
massless in the symmetric phase and heavy in the broken phase), or a tree- level contribution to 
the barrier separating the two phases. In the standard model, the electroweak phase transition is 
not strongly first-order. There are no heavy bosons (relative to the Higgs, which sets the relevant 
scale), and at high temperature the contribution of heavy fermions (top quarks) does not increase 
the barrier since J_(x^) does not contain any cubic terms. 

The particle spectrum in the NMSSM may seem somewhat promising, since there are additional 
heavy masses in the Higgs sector and field-dependent neutralino masses, but these are not enough 
to guarantee a strong transition. Since many more particles couple to the Higgs than to the singlet, 
finite-temperature effects drive {hu) and {hd) to zero at temperatures well below the point at which 
they drive (s) to zero. Therefore, s can be large on either side of electroweak symmetry breaking, 
and some of the new particle masses that depend on s can be heavy even in the symmetric phase. 

However, the NMSSM can succeed in producing a strongly first-order transition through its tree- 
level contributions. If the transition occurs both in the Higgs and singlet directions simultaneously, 
and if the singlet vev is non-zero in the electroweak symmetric phase just above the transition, then 
terms like s'^h'^ and sh"^ both contribute effective cubic terms to the potential which can increase 
the barrier between the the symmetric and broken phases. 



361 ] . We input 



We calculate the phase transition using the software package CosmoTransitions 
the above definition of the effective potential, find the necessary soft-breaking masses that produce 
desired values for tan f3 and fi via a minimization procedure, and choose a renormalization scale 
A such that the one-loop minimum does not drastically differ from its tree-level value. This last 
point requires a certain amount of finesse since the top-quark contribution to the zero-temperature 
one-loop potential tends to be fairly large. The CosmoTransitions package traces the broken 
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A 


0.75 


TO^i [GeV] 


261.26 




0.45 


TO^o [GeV] 


130.72 


tan (3 


1.7 


(crw)f,b [cm^/s] 


3.07 X 10-^^ 


Ax [GeV] 


545.0 


{crv)-y^ [cm^/s] 


1.54 X 10^2^ 


[GeV] 


-88. 


af [pb] 


2.8 X 10-9 


^l [GeV] 


275.8 


flSD [pb] 


1.4 X 10-6 


Ml [GeV] 


143.5 


EWPT Properties: 




Ma [GeV] 


635.5 


Te [GeV] 


72.3 


mft, [GeV] 


126.4 




1.14 



TABLE I. Benchmark Point in the NMSSM with a strongly first- order EWPT and a 130 GeV line. We use 
a renormalization scale of A — 100 GeV in the effective potential. 

electroweak phase up in temperature until it disappears, and then traces the symmetric phase 
down and checks for an overlap. If there is one, it calculates the temperature of degeneracy (the 
critical temperature) and the separation between the phases. If there is no overlap, then the 
transition is necessarily second-order. 

The region of the NMSSM consistent with the 130 GeV Fermi line can in fact accommodate a 
strongly first-order phase transition. The barrier has large tree-level contributions and in particular 
does not require an additional light scalar. As a proof of principle, we outline a benchmark point 
consistent with a 125 GeV Higgs, 130 GeV Fermi line, and a strongly first-order electroweak 
phase transition in Table HI This point has an EWPT at Tc = 72.3 GeV with order parameter 
ip{Tc)/Tc = 1.14 and is consistent with all other relevant phenomenological constraints [f]. The spin- 
dependent and -independent neutralino-proton scattering cross-section for the point in Table H] is 
computed taking a^^f^f = 39 MeV, cjo = 43 MeV for the strange quark content of the proton and 
is thus rather optimistic. Also, note that we do not show the relic density for the specified point. 
Since we are near a resonance, as discussed in Sec. IIIBl the relic density calculation should be 
performed to loop level - something which is not implemented in Mi cr Omegas Neglecting one- 
loop processes, the relic density for this point may be too large. We have checked, however, that 
at tree-level and neglecting the contribution from the resonances, one can introduce a light slepton 
with ~ 140 GeV which will set Qh'^ = 0.11 for the parameters shown. Since tan/3 is small, the 



As mentioned previously, we can invoke some higher-dimension operators to render A perturbative below the GUT 
scale. 

We have also found a suspected numerical issue with the MicrOmegas 2.4.5 calculation of the relic density near 
the resonance. There is a very sharp increase in the annihilation cross section right above tuai ~ 2m^o which we 
believe is unphysical. Since the zero-temperature total-annihilation cross section is of order (av) ~ 10^^® cm'^/s, 



by the arguments in Sec. lIIBTl the thermally averaged cross-section at freeze-out should be smaller than this since 
the resonance is effectively shifted. Instead, we find iSrop of four orders of magnitude in the relic density which 
is quite suspect. 




50 100 150 200 250 300 350 400 450 
T (GeV) 



FIG. 4. The phase structure for the benchmark point with first-order phase transitions. The dotted line 
gives the temperature- dependent singlet field values, and the solid line gives the temperature- dependent Higgs 
doublet field values. 

presence of such a light slepton will not affect the properties of the EWPT. Thus, we are confident 
that a proper one-loop calculation of the relic density for the benchmark point in Table |T] will yield 
a relic density compatible with observation, albeit with some possible changes to the parameters 
or the introduction of a co- annihilation channel which will not substantially affect the EWPT. 

Fig. S] shows the field evolution as a function of temperature for the benchmark point in Ta- 
ble m This makes the location of the phase transitions obvious: first-order phase transitions can 
happen anywhere there is a discontinuous jump in the vacuum expectation values. A second-order 
transition, if there were one, would be distinguished by a continuous line of vacuum expectation 
values with discontinuous first derivatives. 

Fig. [5] shows the field configuration at the critical temperature of electroweak symmetry break- 
ing. All three fields — s, hu and — change values when tunneling from the high-temperature to 
the low-temperature minimum. We calculate the tunneling direction (denoted by a thick black line) 
using the CosmoTransitions package, where by "tunneling direction" we mean the path through 
field space that one would travel when crossing a bubble wall. The path is curved in the s — hu 
and s — hd planes, but is approximately straight in the — plane (A/3 ^ 1). 

While we did not perform a systematic stud y of the NMSSM parameter space compatible with 



a strongly first-order transition (see e.g. Refs. 



IC 



111 ] for previous work in this direction), there 



are some common traits between the viable points we have found. Restricting ourselves to the case 
of positive A, k, /x, and we find that a strongly first-order phase transition typically requires 
A > 0.6, K < 0.6, ^4^ > 500 GeV, and /x < 350 GeV. This seems to be consistent with our intuition: 
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20 40 60 80 100 

vWW (GeV) 

FIG. 5. A contour plot of the effective potential just below the critical temperature. The electroweak broken 
minimum is represented by the dot on the upper-right, while the symmetric minimum is on the lower left. 
The actual tunneling happens along the curved solid black line. 

increasing the strength of the cubic terms in the effective potential and decreasing the singlet vev 
tends to strengthen the transition. Note that, for all the points we considered, the transition tends 
to happen in two steps: the system transitions away from (s) = at a high temperature, around 
300-400 GeV; while electroweak symmetry breaking happens much later, at a temperature around 
or below 100 GeV. 



IV. COMPUTING THE BARYON ASYMMETRY 



The discussion in the previous section makes it clear that a strongly first-order EWPT can 
occur in the NMSSM region of parameter space compatible with the Fermi 130 GeV line. We 
now turn our attention to the CP-violating sources also required for electroweak baryogenesis, and 
to the detailed requirement of producing the correct amount of baryon asymmetry in the early 
universe, parametrized by the baryon-to-entropy raticF^. Yb ~ 10^^*^. As we show in this section, 
CP-violating higgsino-gaugino sources can be very efficient in the NMSSM regions of interest and 
potentially source the observed baryon asymmetry of the universe. 

In electroweak baryogenesis, the baryon asymmetry is produced by SU (2) sphalerons acting on 
a net left-handed chiral density, ni. To determine ul, we must solve a set of quantum transport 
equations for each of the relevant particle densities contributing to the LH charge density. For each 

For concreteness and consistency with previous studies, we take Yb ~ 9.1 x 10~^^ 
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of these charge densities, n,, the Schwinger-Dyson equations yield the continuity equations [37| 

^+V-j-(x) = 5,(x). (21) 

The RHS of the above equation contains both CP-conserving and CP-violating contributions. For 
the case of Dirac fermions, the sources are given by 



S^{x)= d?z dzoTr[S>(x,z)G<(z,x)-G>(x,z)S<(z,2;) 

J J-oo (22) 

+G<{x,z)J:>{z,x) - ^<{x,z)G>{z,x)] 

where G"^'-^ , S^'^ are Green's functions and self-energies, respectively, in the closed time path 
formalism (see e.g. Ref. [sT^ for details). We focus here on the case of gaugino-higgsino sources, 
and compute the quantities Sgo^+in the Higgs vev-insertion approximation, which we describe in 



more detail below (see e.g. Ref. 38|] for a recent discussion on scalar sources in the MSSM). 



A. The VEV-Insertion Approximation 

The CP- violating interactions we consider involve the scattering of higgsinos and gauginos with 
the spacetime-dependent Higgs vevs in the bubble wall. In what follows we parallel the derivations 
for the corresponding quantities in the MSSM found in Ref. js^]- We will assume that the necessary 
CP- violating phase (j) is that of the wino soft SUSY-breaking mass M2 (in fact, the relevant phase 
is the relative phase between Mi^2 and fi, however as discussed previously we take fi, Mi to 
be real to avoid large spontaneous CP-violating effects in the computation of the various dark 
matter properties). The part of the NMSSM Lagrangian giving rise to the relevant CP- violating 
interactions is then given, in terms of four-component spinors, by: 



Ant D -^^HO 



Vd{x)PL + e''t'vu{x)PR -^^,-92^^^+ Vdix)PL + e''f'vuix)PR ^^++h.c. (23) 



- w+ 



where Pl,r are the usual projection operators. 

The spinors ^^o.± satisfy Dirac equations with a spacetime- varying mass fi{x). As discussed 
in Sec. IIIIl the profile fj,{x) depends on the detailed properties of the phase transition at each 
point in parameter space. In the region of interest, however, the singlet vev does not change very 
significantly during the EWPT. Consequently, even though the variation of the singlet vev was 
crucial for achieving a strongly first-order phase transition, we ignore its space-dependence heri^ 
and approximate //(x) by its value after the EWPT, /Li(x) ~ /i. Then the mode expansions for the 

The spacetime-dependence of fi can introduce novel sources of CP-violation in the NMSSM; see e.g. Ref. [l^ 
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operators in the Lagrangian Eq. (\23h are the same as in the MSSM case and so the resulting source 
from Eq. (j22p matches that of the MSSM in the vev-insertion approximation: 



X] iidji^iz) + 9jiz,x)]ReTr 
j=A,B 



G^. , {x, z)G<~^ {z, x) - G%^ (x, z)G>~^ {z, x) 



-i [gj{x,z) -5j(z,x)]ImTr 



G^^ {x, z)G<~^ {z, x) - G^^ (x, z)G>~^ (z, x) 



(24) 



where the sum over B is over contributions arising from momentum and mass terms in the 
spectral function, respectively, and where 



9A{x,y) = — [vd{x)vdiy) +Vu{x)vu{y)] 



gB{x,y) = ^ [vdix)e '^Vuiy) + e'^Vu{x)vdiy) 



(25) 
(26) 



The rest of the derivation proceeds as in the MSSM case, i.e. by performing a derivative ex- 
pansion in gA,B{x,z) around z = x. The CP-conserving sources arise from the terms in Eq. (j24p 
symmetric under the interchange of x -f-T- z and so appear at zeroth order in this expansion, while 
the CP-violating sources arise at first-order. In particular, performing the integration for the 
CP-violating contribution yields 



2 

Sf± = ^vixf$ix)M2tism< 



(27) 



where uj'^- — 
H.w 



Ikl -|- M~ ~ (the masses here include thermal contributions, (Je^), Sn^r = 



iTg^ (here the are the thermal widths of the higgsinos and winos in the plasma) , and 

np is the Fermi distribution function. The corresponding expressions for the CP-conserving (and 
neutral higgsino CP-violating) sources can be found in Ref. 37| with the appropriate replacements. 

The CP-violating source in Eq. (j27p exhibits several important properties. The first term of the 
integrand in Eq. ()27p is resonant for M2 ~ as can be appreciated by rewriting the denominator 
as 



w 



H 



k|^ + (/. + A)2 + <5|-i(r^ + r^). 



(28) 



Thus for a given choice of /U the parameter A determines the strength of the resonance, and 
hence the resulting baryon asymmetry. At finite temperature, fJ-{T) will generally be different from 
/i(T = 0), since the singlet vev varies with temperature. This can be thought of as providing 
a finite temperature correction to A; we neglect this effect in calculating the baryon asymmetry 
across the parameter space, as this difference depends sensitively on the finite-temperature effective 
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potential at each point. Note also that the Fermi distribution functions in the numerator result in a 
suppression of the baryon asymmetry for masses much larger than the electroweak phase transition 
temperature. As an optimistic estimate, we take = 140 GeV in calculating the BAU across the 
parameter space; the SU{2) sphaleron rate (and hence the overall baryon asymmetry) decreases for 
lower temperatures. For example, taking Tc = 100 GeV will decrease the overall baryon asymmetry 
by a factor of about 0.7 across the parameter space (i.e. the CP-violating phase sin(/> at each point 
would increase by a factor of about 1.4). We encourage the Reader to bear this in mind while 
interpreting our results. 

Other important quantities determining the strength of the CP-violating source are the bubble 
wall width (Lw), velocity (vw), and the variation of Higgs vevs across the wall (A/3). This can 
be seen by approximating the bubble wall profile by a step- function, whence /3 w A(3vw/Lw For 
the wall width and velocity we choose Ly^ = 10/T and Vw = .05. These values are inspired by the 
MSSM and will vary depending on the details of the potential and the spectrum for each point in 
parameter space as discussed in Sec. IIIII In our current set-up, since there is only a small degree of 
mixing between Ai and A2, the quantity A/3 to a good approximation scales as in the MSSM, i.e. 
roughly A/3 oc l/m^^ (in our calculation of A/3 we use the full two-loop results of Ref. jssl). Since 
171^2 will vary across the parameter space, A/3 will have an important effect on the parameter space 
available for EWB. For the values of we consider, A/3 falls in the range A/3 ~ 10"^ - 10-"^. 

The other relevant particle number-changing processes (including the triscalar, Yukawa, and 
CP-conserving relaxation interactions) are also coniputed in the vev-insertion approximation; ex- 



pressions for these rates can be found in Refs. 



37|, 



4014431]. In addition to these MSSM processes. 



there are new interactions in the NMSSM arising from the singlet and singlino degrees of freedom. 
In particular, there is a resonant relaxation term (and possible CP-violating source {l^) arising 
from higgsino-singlino interactions with the Higgs vevs. The relevant part of the Lagrangian is 



^int — ^ 



Vu{x)H^S + Vdix)H^S + h.c. (29) 



where ^ and S correspond to the two-component higgsino and singlino fields. We can rewrite 
Eq. (|29|) in terms of four-component spinors as 

Ci, = X^fio[vu{x)PL-Vd{x)PR]^^ + h.c. (30) 



and follow the methods of Ref. [37[ to compute the source. Since we assume that there is no CP- 
violation in the singlino sector, Eq. ()30p results in a resonant chiral relaxation rate for the higgsino 

23 



chemical potential T~n~ = F — - + T~ - where 

^ H^S HOS H^S 



1 A2 /-^ dfcA;2 

— im < t^i ~ — k" — IVl^ sm 2, 



F — E* 

^S 



+ [SgSg + k'^ + Mg\fi\ sin 2/3] 



^5 + 



(31) 



and where the various quantities are defined analogously to those in Eq. (j27p . The singlino mass 
given by 

M| = 4kVVA^ + 4 (32) 
(here 5^ is the singlino thermal mass), and the quantity hp is defined as 

Since we consider moderate values of A, we take ~ O.OOIT for the singlino width. The denomi- 
nator of the first term in Eq. ()3ip has the same resonant structure as in Eq. (I28p and is the most 
significant contribution to the transport equations from the singlino, tending to reduce the resulting 
charge density. Given our choices for A and k in Fig. [3l the relaxation rate F^og is near resonance 
in this region since ~ /i. We account for this higgsino-singlino resonant relaxation in our 
computation of the baryon asymmetry, but do not consider the other non-resonant singlet/singlino 
interactions, as they are subdominant. 



B. Solving the Transport Equations 

With the sources contributing to the RHS of Eq. ()2ip for the various charged current densities 
in place, we compute the baryon asymmetry point-by-point across the 130 GeV line parameter 
space described in Sec. Ill Al for A = 0.6, k = 0.32, and tan/3 = 1.8 as an example. We do so 
by solving the system of transport equations to determine the LH charge density ul, assuming a 
strongly first-order EWPT and that the SU (2) sphaleron rate F^^ is slow compared to the other 
particle number-changing rates. Then, given ni^z), the baryon number density results from the 
integral of hl over the unbroken phase, 

_ — ^£ws_ I nL{z)e ^, (34) 

where z is the comoving distance away from the bubble wall (neglecting the curvature of the wall 
and taking z < to be the symmetric phase). 
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To determine ni, we work under the set of assumptions detailed in Refs. 37|, l4J], and in 
particular assuming "super-equilibrium", i.e. that the chemical potentials of all SM species and 
their superpartners are equal 



42l |. This allows us to define common charge densities for Higgses 
and higgsinos, quarks and squarks, etc. Given the condition of super-equilibrium and that the 
sfermion masses are heavy, one can show that the relevant charge densities we must keep track 
of are those corresponding to the Higgs/higgsinos (H), the right-handed tops/stops (T), and the 
left-handed third-generation quarks/squarks (Q). The transport equations then read 
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(35) 



(36) 



(37) 



Here, Tmt,h are chiral relaxation rates (including the contribution from the higgsino-singlino-vev 
interaction), active only in the bubble walF^ 



tion rates 



and broken EW phase, Tyt are Yukawa interao 



431 ]. Tss is the SU{3) sphaleron rate (responsible for generating densities of first- and 



second-generation quarks), and the fcjS are statistical factors relating the charge densities rij to 
the corresponding chemical ^potential /ij. We solve Eqs. ()35p - (l37p utilizing the diffusion approxi- 
mation discussed in Ref. [37]. The LH charge density entering into Eq. (134p is then given to good 
approximation by the relation 



ULiz) = 5Q{z)+4:T{z). 



(38) 



Contours corresponding to the observed value of the baryon-to-entropy ratio are shown across the 
130 GeV line parameter space on the resonance (A = 0) for different values of the CP-violating 
phase (p in Fig- El 

In interpreting our results, the reader should bear in mind that there are several uncertainties 
present in our calculation of the baryon asymmetry. As mentioned, the microphysical properties 
of the EW bubble wall and details of the electroweak phase transition {L^, f^, A/3, Tc, etc) can 
significantly affect the calculation of and Yb (see e.g. Ref. 4j] and references therein for a 
more detailed discussion of these effects). Also, there are several other frameworks for calculating 



For simplicity, in solving the transport equations we assume a step- function profile for the Higgs vevs in the bubble 
wall. 
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the baryon asymmetry |45l450l| , with results that can differ by up to an order of magnitude from 
one another (for a review of these different approaches, see Ref. [9|). Additionally, there are other 
possible sources of CP-violation in the NMSSM that could contribute to the BAU in this scenario. 
For example, allowing a relative phase between A and k would allow resonant CP-violating singlino 
sources arising from Eq. (pOj) which in fact would be close to resonant (see Ref.p^] for a discussion 
of singhno-driven EWB in the NMSSM). 

Despite these issues and caveats. Fig. [3] suggests that resonant CP-violating higgsino-gaugino 
sources can be very efficient in the region of the NMSSM consistent with a 130 GeV gamma-ray line. 
Even if we had over-estimated the baryon asymmetry by an order of magnitude, there could still 
be regions consistent with both the Fermi line, the observed BAU, constraints from electric dipole 
moments (which we discuss below), and DM direct detection, provided more optimistic choices for 
the strange quark content of the proton or the local distribution of dark matter. For example, 
taking the values of do, cTttN we considered for the EWPT benchmark point pushes out the allowed 
values of Mi in Fig. [3] out to about 145 GeV, which would allow a factor of ten over-estimation of 
the BAU consistent with EDM constraints. 



C. EDM Constraints 

The NMSSM contains several possible sources of CP-violation beyond those in the MSSM: CP- 
violation in tree- level parameters A, k, and fi; CP- violation in soft-breaking terms Ax and A^; 
and additional effects coming from the mixing between the two CP-odd eigenstates Ai and ^2- 
However, in our setup we assume no CP-violation in the tree-level Higgs sector and very little 
mixing between Ai and A2 {Ai must be mostly singlet-like, as explained above). Therefore, the 
electric dipole moment calculations reduce to those in the MSSM. 

We use the package CPSuperH 



5l| to calculate the electric dipole moments of the electron, the 
neutron, and the mercury atom, which have current experimental limits of \de\ < 1.05 x_10~^'^e 
cm 



53|, and \dHg\ < 3 x lO^^^e cm [5j]. The 



5^ (via the YbF molecule), < 2.9 x lO^^Gg 
neutron and the Mercury atom generally provide extremely strong limits on CP-violating physics, 
but they are most sensitive to chromo-EDMs and CP-violation involving colored particles. We have 
no chromo-EDMs in this model, so the electron EDM provides, here, the strongest constraint. All 
one-loop EDMs are suppressed by the heavy sfermion masses. The dominant two-loop contribution 
comes from the Barr-Zee diagram containing a chargino loop. 

For each point in the parameter space of Fig. [3l we calculate the EDMs using the value of (p 
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that produces the proper baryon abundance. Except for most of the parameters necessary for 
calculating the EDMs vary little over the plotted region, so the EDMs are most sensitive to and 
the corresponding iso-level curves follow similar trajectories. The small region in the upper-left 
with sini;^ > 0.37 has \de\ > 1.05 x lO^^^e cm, and is thus ruled out by experiment. The smallest 
EDM in this region, corresponding to smcf) w ^, is \de\ = 5.1 x 10^^®. This is well within the 
anticipated sensitivity of next-generation EDM experiments (for a review, see, e.g., Ref. jssi). 
which have the potential to either rule out or lend credence to this baryogenesis scenario. 



DISCUSSION AND CONCLUSIONS 



The present study reaffirms that the NMSSM framework (and indeed other singlet-extensions 



of the Higgs sector [56(]) can provide a viable explanation of the 130 GeV Fermi gamma-ray line 
in terms of resonant neutralino annihilation through a pseudoscalar into photons. Agreement 
with observation and with the relevant constraints is realized in the NMSSM for a bino-like LSP 
(dictating that Mi ~ 130 GeV), with relatively large A, moderate /i, and with Ai predominantly 
singlet-like to avoid indirect detection constraints on continuum photons. While there are many 
independent constraints on this scenario, currently there remains a substantial amount of parameter 
space consistent with the gamma-ray line and in agreement with the various dark matter and 
particle physics constraints. 

Here we have shown that the parameter space consistent with the Fermi line in the NMSSM 
is also promising for electroweak baryogenesis. In particular, the relatively large values of A typi- 
cally considered tend to bolster the cubic term in the finite-temperature effective potential in the 
direction of electroweak symmetry breaking, leading to a strongly first-order electroweak phase 
transition in parts of the parameter space. Additionally, the moderate values of /U ensure that the 
singlet vev is not too far from the EW scale, again tending towards a strongly first-order transi- 
tion. We illustrated this in Sec. IIIII by providing a benchmark point consistent with the 130 GeV 
line and a strongly first-order EWPT, and in agreement with all other relevant phenomenological 
constraints. While we only studied in detail one particular point as a proof of principle, we expect 
a more systematic study of the NMSSM parameter space to uncover many other regions consistent 
with the line and a strongly first-order EWPT. 

Not only does the parameter space consistent with the line support the possibility of a strongly 
first-order transition, it can also provide an efficient source for CP-violation that gives rise to the 
observed baryon-to-entropy ratio of the universe. Resonant higgsino-gaugino sources can be very 
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efficient here due to the moderate values of Mi^2 and ^ required to produce the Une. In particular, 
allowing for a CP-violating phase in M2 does not strongly affect the line or the dark matter 
phenomenology, but it can produce the observed BAU with sin (f) small enough to be consistent with 
electric dipole measurements, as shown in Sec. IIVI and Fig. [3l While we focused on the higgsino- 
wino sources in the present study for the sake of illustration, similar resonant CP-violating sources 
arising from other interactions can be active in the same regions of parameter space by similar 
reasoning. For example, if one allows for Mi to carry a complex phase, resonant bino-higgsino 
sources can be very efficient as well. This may be of particular interest in the case of negative 
fi whereby |^| can be taken as low as 140 ~ 150 GeV (and thus potentially very close to this 
resonance) while in agreement with direct detection constraints jl?! . A careful study of the effect 
of a CP-violating phase in Mi on the line and dark matter properties would be necessary to assess 
whether such a scenario is possible, but we expect it is since EDM measurements dictate that the 
CP-violating phase is necessarily small. Also, singlino-higgsino sources can in principle be efficient 
in this region as well, provided a relative phase between A and k [12], again due to the moderate 
values of the singlino mass (see Eq. (j32p ) and fi in this scenario. These other sources would be 
especially important for points such as our EWPT benchmark which features a rather heavy wino 
but lighter bino and singlinql^. 

An interesting feature of our scenario is that the relevant parameter space will be conclusively 
tested in the near future by modest improvements in various experimental efforts. The moderate 
values of fi we consider result in rather large spin-independent neutralino-nucleon cross-sections 
which continue to be probed by direct detection experiments. The relatively large values of A, as 
required for a large (cru)^^, combined with the large A\ and moderate values of k necessary for 
a strongly first-order EWPT, tend towards a significant coupling of Ai to e.g. bb and so will be 
tested by modest improvements in indirect detection experiments. Additionally, the CP-violating 
phase(s), required to source the left-handed charge density for the SU{2) sphalerons, will be well 
within reach of various future EDM experiments (see e.g. Ref. 3a] for a related discussion). The 
whole scenario will also continue to be tested by ongoing measurements of Higgs couplings and 
searches for other particles at the LHC. 

Of course the viability of the 130 GeV line scenario in the NMSSM or any other SM extension 

hinges on the persistence of the line in the Fermi data and on a dark matter interpretation of these 

results. If the line is indeed due to resonant dark matter annihilation, this work shows that the 

^'^ Note that non-resonant wino-higgsino sources, such as those considered in Refs. S) 0i '^^^ ^^^'^ potentiaUy 
provide the necessary CP- violation for our particular EWPT benchmark point. 
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NMSSM framework can potentially explain the origin of both the baryonic and dark matter in our 
universe. 
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Appendix A: Scalar Mass Terms 

We present here the scalar mass terms used in the calculation of the finite-temperature effective 
potential. These are simply the second-derivatives of the full 10-parameter potential, but simplified 
such that only 3 of the parameters and s) are non-zero. Each subscript denotes a partial 

derivative with respect to that field. Primed subscripts are derivatives with respect to the imaginary 
field components, and u and d denote derivatives in the up- and down-type charged directions. The 
tree-level masses are just the mass eigenvalues of the following matricies. 



ML =hKhl + s') + Ugl + giml - hi) + ml (Al) 

Mjd =\Khl + s^) + \{gl+ gl){2>hl - hi) + ml (A2) 

=\X{hl + hi) + 3k^s^ - Xnhuhd + ml + V^nA^s (A3) 

Mid =>?Khd - |Aks2 _ 1(^2 ^ gl)h^hd - -^XAxs (A4) 

Ml, =\^hus - Xnhds - -^XAxhd (A5) 

Ml =\^hds - Xnhus - -^XAxhu (A6) 

Mlu' =\Khl + s^) + \{gl + gl){hl - hi) + ml (A7) 

Mid' =\Khl + s') + Ugf + gl){hl - hi) + ml (A8) 

Myy =^X{hl + hi) + K^s'^ + Xnhuhd + ml - \P2kA^s (A9) 

mI,^,=\Xks^ + ^^XAxs (AlO) 

Ml,, = ~ XKhdS + -j^XAxhd (All) 

Mj,,, = - XKhuS + -^XAxhu (A12) 
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Mlu =h>?s^ + \{gl + gl){hl - hi) + \glhl + ml (A13) 
^i- =^^'^' + + 9l){hl - hi) + \glhl + ml (A14) 
^Id = ~ l^^^uhd + ^Aks^ + \glhuhd + ^A^as (A15) 

There is a second matrix for the charged Higgs, but the two are identical except for a change 
of sign in the off-diagonal term which does not affect its eigenvalues. 

In the high-temperature approximation, the thermal mass terms come from the quadratic piece 
of the one-loop finite-temperature contributions to the effective potential. The scalar thermal 
masses include contributions from all particles with field dependent masses, and they get added to 
each of the diagonal terms in the mass matrix. They are: 



n, = r2 [\ {gl + 2,gl) + \X^ + \yl] 

n, = ir2(A2 + K2). 



(A16) 
(A17) 
(A18) 



The longitudinal polarizations of the gauge bosons also receive thermal mass corrections. At 
finite temperature, the gauge boson mass mixing is 

^ {u ^ 

hl + hl 



M, 



gauge— long 



9l 



V 



92 9192 

^2 



hi 



\ 



(A19) 



6 9i/ 



9192 9i j 

Again, we have ignored the contributions from the sfermions, because they are much too heavy 
to factor into the high-temperature corrections. For more information on thermal masses in the 
supersymmetric theories, see ref. [58i] . 
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